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ABSTRACT 

Whistler propagation in the magnetosphere is studied in detail to 
find accurate and economical means of determining the path latitude and 
the electron concentration along the path from whistler parameters f 

n 

(nose frequency) and t (travel time at the nose) . Longitudinal propaga- 

n 

tion in field aligned whistler ducts of cold plasma is assumed, and the 
earth’s magnetic field is approximated by a centered dipole. The effects 
of whistler propagation in the earth-ionosphere waveguide and through the 
conjugate ionospheres are treated as small perturbations. Several alter- 
nati ye methods are described so that the most economical method may be 
chosen depending on the desired accuracy and the availability of a com- 
puter or a calculator. 
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I . INTRODUCTION 


A. INTRODUCTION 

When Storey successfully explained whistlers in 1953 as radiation 

from lightning propagating from hemisphere to hemisphere along geomagnetic 

lines of force, he also discovered that an equatorial electron concentra- 
-3 

tion of ~ 400 cm at a geocentric distance of 3 R was required to explain 

E 

the observed whistler dispersion [Storey, 1953]. This was the first evi- 
dence of such large electron concentrations far beyond normal F2 layer 
heights, and it opened a new region of interest for exploration which 
later came to be called the protonosphere, magnetosphere or plasmasphere . 
Since then whistlers have provided a powerful and unique tool for probing 
this remote region of space, and much progress in our understanding of the 
morphology and dynamics of the inner magnetosphere has been made through 
grour.d-based whistler observations. 

The discovery of nose whistlers [Helliwell et al . , 1956] added the 
possibility of determining whistler path latitude and permitted more de- 
tailed studies of plasma distribution in the magnetosphere [Smith, I960], 
This led to the discovery of the plasmapause [Carpenter, 1962, 1963] and 

later to detailed descriptions of its dynamical behavior and to measure- 
ments! of electron concentration profiles inside and outside the plasma- 
pause [Carpenter, 1966; Angerami and Carpenter, 1966; Carpenter, 1970; 

Park and Carpenter, 1970]. More recently, the whistler technique has 
been used to measure cross-L motions of thermal plasma in the plasmasphere 
[Carpenter and Stone, 1967; Park and Carpenter, 1970; Carpenter et al,197l] 
and coupling fluxes of ionization between the protonosphere and the iono- 
sphere [Park, 1970]. 
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As the whistler technique finds more applications in magnetospher ic 
research, there is a growing need for simple means of processing the data 
and for improvements in accuracy. This report describes various ways of 
deducing information on path latitude and electron concentration from 
whistler nose frequency and travel time at the nose. The whistlers are 
assumed to be received on the ground. When the nose is not visible on a 
spectrogram, the nose frequency and travel time at the nose can be extra- 
polated from the visible portion of the whistler trace by the techniques 
such as those developed by Smith and Carpenter [1961], Dowden and Allcock 
[1971] and by Bernard [1971]. 

This work is an extension of the work of Angerami [1966], who refined 
the whistler technique by obtaining empirical support for a diffusive- 
equilibrium model of the field-line distribution of electrons inside the 
plasmapause and an idealized collisionless model outside the plasmapause. 

The basic approaches and formulations used here are similar to his, and 
many of the ideas for improvements developed out of conversations with him. 

For the theory of whistler propagation and other background material, 
the reader is referred to a book by Helliwell [1965] or a review paper by 
Carpenter and Smith [1964] . For a detailed development of models of 
electron distribution along geomagnetic field lines, see Angerami and 
Thomas [1963] and Angerami [1966]. 

Figure 1 shows an example of a multi-component whistler train. At 
the top is a frequency-time spectrogram of a recording from Eights , Antarctica . 
The causative sferic is identified and marked t = 0. The middle panel 
shows a tracing of three nose whistler components, while the bottom panel 
is a sketch of their field line paths. The purpose of this report is to 
develop methods of determining accurately and economically the path latitude 
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Figure 1. An example of a multicomponent whistler recorded at Eights, 
Antarctica on June 6, 1963. The top panel is a frequency- 
time spectrogram showing many discrete whistler components all excited 
by a sferic at t = o. The middle panel is a tracing of 3 whistler com- 
ponents and the causative sferic. The bottom panel is a sketch of cor- 
responding field-aligned propagation paths in the magnetosphere. 
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and electron concentration along the path from whistler parameters f and 

n 

t . 

n 

B. THE WHISTLER INTEGRAL 

In the case of longitudinal propagation of a whistler wave, the refrac- 
tive index p is given by 

1/2 



where f is the wave frequency, f the plasma f requency , and f^ the 
electron gyrof requency . If p » 1, calculations can be greatly simpli- 
fied by neglecting the 1 inside the bracket. This is a good approximation 
inside the plasmapause, where plasma frequencies are relatively high. For 

_3 

example, assume an equatorial electron concentration of 100 cm at L = 4 

(f = 90 kHz and f =13.5 kHz.) At f = 0.5 f , p ^ 25 and the error 
p H H’ 

introduced by neglecting the 1 is less than 0.1%. In the plasma trough 

region beyond the plasmapause, the error is larger because of low plasma 

densities there. Angerami [1966] examined in detail the validity of this 

approximation when applied to ’knee whistlers’ propagating just outside 

the plasmapause. He found that on an occasion when the equatorial electron 

-3 

concentration was as low as 2 cm at L = 4, the above approximation 
introduced an error of ~ 40% in the calculated electron concentration. 

When electron concentrations are that low, errors due to other factors 
which are negligible under normal conditions also become serious (see 
Chapter 6). For extremely low electron concentrations, all sources of 
error including the 1 in Eqn. (1.1) should be carefully evaluated. 

The refractive index is then written as 
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\1 = 



1 


and the group refractive index as 


= df ( " f) = 


2f 


V 2 f l/ 2 /!_ 1_\ 3 / 2 


H 




The "ravel time of a whistler wave propagating along a geomagnetic field 


line is 


t(f) 


1 _ 

2c 




T[2 


ds 


( 1 . 2 ) 


where c is the speed of light in free space and ds is an element of 

path length along the field line. The integral is over the entire length 

of the field line path. If a centered dipole approximation is used for 

the =>arth T s magnetic field, f is a simple function of s (see Appendix) . 

H 

In order to express f^ as a function of s, however, it is necessary' to 
adopt a model of the plasma distribution along field lines. Since the 
exact distribution along a field line at a given time is not known, the 
choice of a model is somewhat arbitrary. The results, however, are re- 
markably insensitive to the choice of models as will be seen in the next 
chapter . 

The solutions of Eq . (1.2) show that t has a minimum value t at 

n 

some frequency f for all reasonable models of f . In other words, Eq . 
n p 

(1.2) describes the shape of a nose whistler. Our task is to establish 

quantitative relationships between (f , t ) and (f , f ) , so that the 

H n n H p 

infcrination in the whistler path latitude and on electron concentration 
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along the path can be determined from observed whistler parameters. The 

integral in Eq . (1.2) can be evaluated analytically only for extremely 

l/2 

simple models of f , i.e. f = constant or f oc f ' (gyrof requency 
p p pH J 

model). For more realistic models adopted in this report, numerical 
techniques must be used. 

Because of the gyrof requency term in Eq . (1.2), electrons at low 
altitudes such as in the ionospheric F2 region make a relatively small 
contribution to the total integral, in spite of their large concentrations. 
It will be shown in Chapter 6 that if the plasma is in diffusive equili- 
brium along magnetic field lines, nearly 80% of the propagation delay 
occurs within 30° of the magnetic equator. For this reason, whistlers 
can be used to measure electron concentrations near the equatorial plane 
with great accuracy, without a precise knowledge of concentration profiles 
in the F region and the topside ionosphere up to several thousand kilo- 
meters. Further discussion on this point follows in Chapter 6. 

C. AN OUTLINE OF STRATEGY 

This section contains an outline of the methods of attacking the 
problem and of the organization of the remainder of this report. 

As pointed out in the previous section, the ionospheric F2 region, 
with its large electron concentrations and great variability, makes a 
relatively small contribution to the whistler integral. Therefore, the 
conjugate path segments through the F region will be treated as a small 
correction term. An arbitrary boundary is drawn at 1000 km altitude, and 
the path of integration in Eq . (1.2) is divided into the T ionospheric T 
part below 1000 km and the * magnetospheric f part above (this will serve 
as the definition of the ionosphere and the magnetosphere throughout this 
report, unless otherwise stated). Equation (1.2) is then rewritten as 
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t(f) 


1 _ 

2c 


ds 


/ 


mag 


V 2 f . V 2 /!- *_\ 3 / 2 

H ) 


H 


ds 




lono 




372 


(1.3) 

The two integrals are illustrated schematically in Figure 2. The dashed 
curves 1 and 2 represent the ionospheric and the magnetospher ic contri- 
butions, respectively, and the solid curve 3 is the sum of the two. The 
relative contribution by the ionosphere is exaggerated for the purpose of 

illustration. An actual whistler with observed f and t would have f’ 

n n n 

and r f if it had traveled through the magnet ospheric path only, 
n 

In Chapter 2, the ionosphere is ignored, and methods are developed 

to obtain the whistler path latitude and electron concentration along the 

path from f^ and t^. In Chapter 3, the ionospheric part of Eq . (1.3) is 

considered, and methods are developed to obtain f' and t* from f and t . 

n n n n 

If t is measured from the sferic on a whistler spectrogram (see Figure 34) , 
n 

the jsubionospheric propagation time must be taken into account. This is 

discussed in Chapter 4. Chapter 5 summarizes the results with step by 

step instructions for calculating magnetospheric parameters from f and 

t . In Chapter 6, uncertainties in the whistler method are discussed in 
n 

detail and a few suggestions are made for future improvements. Finally 
for purposes of easy reference an appendix presents a collection of formulas 
and graphs involving some frequently used parameters of a dipole geomagnetic 
field . 
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II. CALCULATIONS AND RESULTS 


A . OU TLINE OF THE CHAPTER 

In this chapter, we consider whistler wave propagation in the magneto- 
sphere from 1000 km altitude in one hemisphere to the 1000 km level in the 


conjugate hemisphere. We assume that the magnetosphere is symmetrical 


about the equator and write 


equator 

,( f , = i / 

1000 km 



( 2 . 1 ) 


Several different models for the distribution of f along field lines will 

P 

be adopted in the next section. With a model for f , the integral in Eq . 

P 

(2.1) can be evaluated numerically for a given magnetic field line. The 

wave frequency f is varied in search for the frequency f ’ which gives the 

n 

minimum time delay t f . This is repeated for different magnetic field lines 

n 

(L = 2, 3, 4, etc., L being Mcllwain’s parameter [Mcllwain, 1961] and i or 

different models of plasma distribution along the field lines. From these 

results, relationships will be found between the whistler parameters f' 

n 

and t’ and the parameters of the medium, L and f . The results wiLL be 
n p 

presented in various forms so that, depending on the desired accuracy, the 

whistler path latitude and the electron concentration along the path can 

be determined by the most economical means. 

We will briefly examine the properties of Eq . (2.1) and justify in 

advance the form in which the results are presented in Sections 2-D, K 

and G, Calculations in the next section will show that f T is approximately 

n 

proportional to the minimum electron gyrof requency along the propagation 

path, f , so that 
Heq 
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( 2 . 2 ) 


f = Kf f 
Heq n 


where K is a quasi-constant that varies slowly with f 1 and models of 

n 

f . The L-value corresponding to f for a dipole field can be easily 
p Heq 


obtained from the relation, 


( f \ 1/3 
Hoeq \ ' 


8.736x10 

f 


(f is the electron gyrof requency at the earth’s equator and its 

Hoeq 

numerical value given above corresponds to a magnetic field strength of 

0.312 gauss.) Since the bulk of the whistler propagation delay occurs 

in the equatorial region where f and f are slow functions of s, Eq . 

p H 

(2.1) can be roughly approximated by a simple algebraic equation, 


t(f) = 


f peq S 1 

c f v 2 7 f 

Heq [ 1- 

\ Heq 


where f is the plasma frequency at the equator, and S_ the path length 
peq 1 

between 1000 km altitude and the equator. At the nose frequency f T , it 

^ ^ t 

follows from Eq . (2.2) that the quantity inside the parentheses 1 - - 

Heq 

is very nearly constant. Since is roughly proportional to L and 

-3 

f oc L , we can write 
Heq 


t’(f') OC 

n n 


5 \ 1/2 


where n is electron concentration at the equator. For a given field line 
eq 
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distribution model, n is proportional to electron concentration at 1000 
km altitude n^ . Therefore, we can also write 


t'(f * ) 

n n 



l/2 


The above two equations can be written 


and 


eq 


= K 


eq 


f f t T 
n n 


f ! t • 
n n 


(2.3) 


(2.4) 


where K and K are quasi -cons tants varying slowly with f’ and electron 
eq 1 * n 

concentration models. 

A very useful parameter in studying the magnetosphere is tube content, 

the total number of electrons in a tube of force. For the purpose of 

2 

defining tube content the tube of force is chosen to have 1 cm cross- 
sectional area at 1000 km altitude and to extend to the dipole equatorial 
Pi aie. Tube content is then 


N T = f n A ds 


where n is the electron concentration and A the cross-sectional 

the tube of force. The integral extends from 1000 km altitude 

equator. Since most of the tube volume is near the equator, we 

tube content to be roughly proportional to n A S. Since S 

eq eq 

1 3 

proportional to L, and A oc oc L , N varies roughly as 

eq Heq 


area of 
to the 
expee t 
is rough 1\ 


n 

eq 


L . 
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Using Eq . (2.3), we write 


2 


N 


T 


K 


T 


f T t 
n n 

~ 


(2.5) 


where K t is again a quasi-constant. In Eq . (2.2) through (2.5), the K’s 
are nearly constant so that once they are evaluated accurately as a func- 
tion of f’ and for various electron concentration models, their value can 
n 

be read from graphs or approximated by empirical formulae with a high de- 
cree of precision. 

B. ELECTRON CONCENTRATION MODELS 


Four different models of electron distribution along magnetic field 

lines will be used in the next section. Angerami and Carpenter [1966] 

compared whistler results with electron concentrations at 1000 km altitude 

measured by the topside sounder aboard the Alouette I satellite, and found 

that a diffusive equilibrium model is a good approximation inside the 

2 3-3 

plasmapause where equatorial electron concentrations are 10 -10 cm 

Outside the plasmapause, where equatorial electron concentrations are 
0 1 -3 

^ 10 -10 cm , collisions may not be sufficiently frequent to establish 
a diffusive equilibrium distribution. Angerami and Carpenter [1966] showed 
that outside the plasmapause, an idealized collisionless model was more 
appropriate than diffusive equilibrium model. However, the choice of an 
appropriate model for low concentration regions remains to be put on firm 
ground . 

The idealized collisionless model to be used in this report probably 
represents an extreme case as compared to most circumstances that occur 
outside the plasmapause. The actual distribution may fall somewhere be- 
tween a diffusive-equilibrium and collisionless distribution. Because of 
the dynamic nature of the plasmapause region, some flux tubes undergo 
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large fluctuations in electron concentration, and the field line distribu- 


tion must change accordingly. For example, during magnetic storms when 
the plasmapause moves inward, a collisionless model may temporarily be 
adequate to describe the low concentration region outside the storm-time 
plasniapause . During the recovery phase, the depleted flux tubes fill 
gradually from the ionosphere, and the field line distribution must become 
more like a diffusive equilibrium distribution as the electron concentra- 
tion increases. For a distribution that might exist during the transition, 
we w:.ll consider a hybrid model in which the plasma is assumed to be in 
diffusive equilibrium from 1000 km altitude up to 30° in dipole latitude, 

but no collisions are allowed within 30° of the equator. 

-4 

In the fourth model to be considered, n cc r along a lield line, 
where r is the geocentric distance. As pointed out by Angerami [1966], 
this model closely approximates the collisionless model adopted here, and 
it has the advantage of computational simplicity. There now follows a 
description of each model and the abbreviations by which each model will 
be identified throughout this report. 

1. Diffusive Equilibrium (DE model) 

Detailed treatment of the problem of a mul t i -component plasma in 
hydrostatic equilibrium along geomagnetic field lines can be found else- 
where in the literature (see, for example, Angerami and Thomas [1963] or 
Angerami [1966]). Here we adopt Angerami’s [1966] formulation as given by 
his Eqs. (2.4) through (2.7).* The electron concentration is given by 

*Angerami stated in a private communication that his Eq . (2.7b) is in 
error and should read 


a ~ r 

o 


R 2g o \ 


2 2 2 
R cos 0 - cos 9 
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n 


( 2 . 6 ) 



z 




H. 

1 



(2.7) 


( 2 . 8 ) 


where 

r = fractional abundance of ionic species 
r = geocentric distance 

Q = angular rotational speed of the earth 
g = acceleration of gravity 
f) = dipole latitude 
k = Boltzman's constant 
T = temperature 
m = mass 

The subscript i refers to the ith ionic species (0 + , He + and H + ) , and 
the subscript 1 refers to the reference level at 1000 km altitude. Figure 
3 illustrates the coordinates and the symbols used in Eqs . (2.6) through 
(2.8). The geopotential height z takes into account variations of gravi- 
tational force with distance and the centrifugal force due to corotation 
of the plasma with the earth. The centrifugal force term becomes impor- 
tant only for L > 6 [Angerami, 1966]. For simplicity, it is assumed that 

the magnetospheric plasma corotates with the earth at all latitudes of 
interest here. At high latitudes, particularly outside the plasmapause, 
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this may not be a good assumption. The effects on the field-line distri- 
bution of plasma of centrifugal force associated with convective motions 
in the outer magnetosphere need further investigation. 

Equations (2.6) through (2.8) show that electron concentration de- 
pends on the temperature and the ionic composition at the reference level 
as well as on geometric factors. We will adopt several different combina- 
tions of the temperature and ionic composition at 1000 km altitude as 

summarized in Table 1. It is assumed that T = T = T. and that T is con- 

e l 

stant with altitude. Since whistlers essentially measure electron con- 
centrations near the equator, it is preferable to normalize n to the equa- 
torial concentration rather than to the concentration at the base of the 
magnetosphere. Equation (2.6) is rewritten as 


n 


n 

eq 


1 hi exp 

2 exp 
i 1 



1/2 


(2.9) 


where the subscript eq refers to the dipole equator. Figure 4 illustrates 
how electron concentration varies along lines of force for different L- 
values according to the DE-1 model (see Table 1) . The distance along the 
field line S is measured from the equatorial plane, and the terminal point 
of the curves correspond to the base level at 1000 km altitude. Near the 
base level, 0 + and He + concentrations decrease rapidly with altitude and 
H + soon becomes the major ion. Beyond that point, electron concentration 
varies slowly with distance because of the large plasma scale height. If 
the temperature or ionic composition is changed, the curves in Figure 4 
remain essentially unchanged at distances more than a few hundred kilo- 
meters above the base level. The only significant changes in electron 
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TABLE I. PARAMETERS OF DIFFUSIVE EQUILIBRIUM MODELS 





Composition 

at 

1000 km. 

Abbreviation 

Temperature 

4- 

O 

H + 

He + 

DE- 

1 

1600°K 

90 r r 

8°r 

2 f 7 

DE- 

2 

3200 

90 

8 

2 

DE- 

3 

1600 

50 

40 

10 

DE- 

4 

800 

50 

40 

10 
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distribution occur where 0 + and He + concentrations are important. 

2 . Collisionless Model (CL model) 

We adopt the formulation of Angerami [1966] for the distribution of 
a proton-electron plasma along lines of force above 1000 km altitude under 
the influence of gravitational force, a charge separation electric field, 
and centrifugal force due to corotation with the earth, but not the effects 
of collisions between particles . The concentration is given by 



z 


2H(l - 



( 2 . 10 ) 


wher? B is the magnetic field strength and z and H are given by Eqs . (2.7) 
and (2.8), respectively. Eliminating n^ in Eq . (2.10), 


n 


n 

eq 



( 2 . 11 ) 


We? assume a constant temperature of 3200°K for both electrons and pro- 
tons . Figure 5 shows the normalized electron concentration as a function 
of distance? along lines of force for several L-values. The CL model 
differs from DE models in presenting a much faster variation in electron 

concentration with distance, particularly near the equator, 

-4 

3 . n oc r Model (R-4 model) 

This model approximates the collisionless model very closely, and 
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20 


Figure 5. A plot of normalized electron concentration in a format similar to that of 
Figure 4 but according to an idealized collisionless model. 




has tie advantage of being much simpler. 


4 -4 

n - n r r 
1 1 


or 


4 -4 

n = n r r 
eq eq 


CM2) 


(2.13) 


4 . Hybrid Model (HY model) 

This model combines a diffusive equilibrium model and a collision- 
less model. We apply DE-2 model (90% 0 + , 2% He + and 8% H + at 1000 km 

and T = 3200°K) from 1000 km altitude up to 30° dipole latitude and CL 
mode] (collisionless H + and electrons at 3200°K) from 30° to the dipole 
equator. This is an attempt to recognize the fact that when the magneto- 
spheric concentrations are low, the plasma may not be in diffusive equi- 
librium, but at the same time, it is unrealistic to assume no collisions 

down to 1000 km altitude, as is done in the CL model. The electron con- 

centration according to this hybrid model is illustrated in Figure 6. 

The models described above are listed on page xi for quick reference. 
C. CALCULATIONS 

Eq . (2.1) is rewritten as 


t(f) 


l 

c 



f 

P 


wv /2 < i 



3/2 



(3.14) 


where ■ 0 is the magnetic 
alti tude . For a dipole 

2 , 

r = r L cos 0 
o r 


latitude, and 0^ the magnetic latitude at 
field (see Appendix) , 


and 


dr 

d0 


2r L cos$ sine) 
o 


1000 km 
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S (x 10 km) 


Figure 6. A plot of normalized electron concentration 
in a format similar to that of Figure 4 but 
according to a hybrid model (see the text for details) 





so that 


ds 

d0 


/dr\ 2 2 

(d^) + r 


V 2 I 2 \ 1 / 2 

= r^L cos0 (1+3 sin 0] 


(2.15) 


The electron gy rof requency can be expressed as (see Appendix) 


f = f 
H Hoeq 


(l + 3 sin 2 0) 1/2 


(2.16) 


where ^ Hoe q is the electron gy rof requency at the earth’s equator, or 


8.736 x 10* Hz. The plasma frequency is 


where a, 
mod 6? Is 
tion 0 


f = ot 

p 


1/2 

n / 


a 


n „ \ (0) 

eq r 


1/2 


is a constant, and \(0) represents various electron 
described in the previous section. The upper limit 
is obtained from the relation 


(2.17) 

concent rat i on 
in the integra- 


. , r i \ 2 / 2 

° os ' 4 i • — 


(2.18) 


Using Eq. (2.15) through (2.17), Eq . (2.14) can be written as 


K 


<f) = 6 " eq 1 ^ 2 / 3((6. Odft s e %‘ /2 K f ) < 2 -i 9 > 

o 


where p is a constant. 

.'For a given field line, say L = 4, 0 is evaluated to the nearest 
even tenth of a degree (the maximum error in 0^ is then 0.1°) . The in- 
tegral is evaluated numerically by Simpson’s rule with /\0 = 0.1°, and by 
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iteration, we search for the wave frequency f^ which gives the minimum 

value of the integral I(f ! ). The iteration is continued until the error 

n 

in f is reduced to less than 0.2%. This procedure is repeated for various 
n 

L- values and various electron concentration models. The results are tab- 
ulated in the second column of Tables 2 through 8 in the next section. 
Following the reasoning given in Section 2A, we express the numerical 
results in the same form as in Eqs . (2.2) through (2.5). The values of 
various K T s are also tabulated in Tables 2 through 8. 


Eq. (2.19), 


eq 


BKf;) 


when expressed in the form of Eq . (2.3), yields 


K 


eq 


3Kf’) 


The electron concentration at 1000 km altitude is easily calculated from 

n and the electron concentration model \(0) : 
eq 


f T t ! 

”l = n eq '■'V = K 1 ~~T 

L 


K and K are then related by 
1 eq 

K 1 = K eq 

The tube content is 


K 


N = 
T 




A ds ,a 
n A tt d0 


d0 


( 2 . 20 ) 


where A is the cross-sectional area of a tube of force with 1 cm area at 
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TABLE 2. RESULTS FOR THE DE-1 MODEL 
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TABLE 3. RESULTS FOR THE DE-2 MODEL 
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TABLE 4. RESULTS FOR THE DE-3 MODEL 
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TABLE 5. RESULTS FOR THE DE-4 MODEL 
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TABLE 6. RESULTS FOR THE CL MODEL 
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TABLE 7. RESULTS FOR THE R-4 MODEL 
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TABLE 8. RESULTS FOR THE HY MODEL 
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1000 km. It is inversely proportional to the magnetic field strength and 


can be written as 


B 1 

JL .\ 3 

‘ 1 

+ 3 

• 

sm p^ 

1/2 

cos0 

6 

1 

2 , n 

+ 3 sin 0^ 

B 


1 

+ 3 

sm p 


cos0^ 


r 

+ 3 s i n0 


is given by Eq . (2.15). 
Eq. (2.20), 


Substituting 


these and n = n 

eq 


A ( 0 ) 


into 


N = 


n r L 
eq o 


(l + 3 sin 2 0 j 


1/2 




~xr 

"OS <p^ 


J* \(0)cos 7 0d0 (2.21) 


The integral is evaluated numerically by Simpson’s rule with \$ = 0.1°, 
and the results are expressed in the form of Eq . ( 2 . 5) , 


N m = 
T 


K 


f ’t’ 
n n 


(2.5) 


D. TABULATION OF RESULTS 

The results of the calculations are listed in Tables 2 through 8 for 

various electron concentration models. These tables represent numerical 

relationships between (f ’ , t T ) and (L, n , n , N ) for selected L-values. 

n n eq 1 T 

The tables can be stored in a computer, and corresponding relationships 
for intermediate values of L can be calculated by interpolation. Para- 
bolic interpolations involving t^ and log (JM) are usually satisfactory. 

E. GRAPHICAL PRESENTATION OF RESULTS 

In this section, the results of calculations are presented in various 
graphical forms whose usefulness depends on desired speed and accuracy. 

In Figure 7, is plotted versus L for a diffusive equilibrium model and 
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Figure 7. A plot of versus L and invariant latitude J\ for diffusive equilibrium and 
collisionless models of field-line plasma distribution. 



collisionless model. The differences among various DE models and the 
difference between the CL model and the R-4 model cannot be resolved on 

plots of this scale. The curve for the hybrid model would lie about mid- 

way between the two curves shown. 

Figures 8 through 11 can be used for a quick determination of n 

eq 

and from f^ and t^. Figures 8 and 9 represent a diffusive equilibrium 
model, and Figures 10 and 11 a collisionless model. Again, the differences 
among various DE models and the differences between the CL Model and R-4 

Model are too small to be shown in plots of this scale. It is not prac- 

tical to obtain n directly from f’ and t* in the manner of n and N . 

1 n n eq T 

The reason is that n^ depends sensitively on the electron concentration 

model, and thus will require a large number of graphs, each corresponding 

to a different model. An alternative approach is to obtain n from n 

1 eq 

or N^. Figures 12 and 13 show the relationships between n , n^ and 

for various electron concentration models. 

As pointed out in Section 2A, a greater precision can be achieved by 

the use of quasi-constants K, K , and . Figures 14 through 21 show 

plots of various K T s versus f’ for various electron concentration models. 

n 

These values of K's are then used in Eqns . (2.2) through (2.5) to calcu- 
late L, n , n 1 , and N T . 

F. OVERLAY METHOD 


For a quick estimate of L, n and N directly from a whistler spec- 

eq T 

trogram, transparent overlays have been devised as illustrated in Figures 
22 and 23. These overlays are for a diffusive equilibrium model, but 
similar overlays can be easily constructed for other models from the re- 
sults in the previous section. The differences among various DE models 
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10"' 10° I0 1 I0 2 

DE MODEL 


Figure 8. Curves of constant equatorial electron concentration for a diffusive 

equilibrium model plotted in t' - f 1 coordinates. The scale at top 
. n n 

shows L-values. 
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CL MODEL 


Figure? 10. Curves of constant equatorial electron concentration for a collision- 
less model. The format is similar to that of Figure 8. 
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CL MODEL 


Figure 11. Curves of constant tube content for a collisionless model. 
The format is similar to that of Figure 8. 
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Figure 13. A plot of the ratio between the tube content and 

the equatorial electron concentration as a function 
of L for different field-line distribution models. 
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Figure 16. A plot of K versus f 1 . DE-1 - DE- 











































Figure 19. A plot of K versus 
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Figure 21. A plot of K versus 
















represent constant values of tube content assuming diffusive 
equilibrium distribution of plasma along field lines. 



canno; be resolved on figures of this scale. The horizontal bars at the 
edges of the figure show L-values, and the sloping curves represent con- 
stant values of n and . These overlays are designed to be used on 
35 mm Rayspan records with 0-20 kHz per 30 mm vertical scale and 20 mm 
per second of real time horizontal scale. An overlay may be placed over 
a spectrogram with the t = 0 line on the causative sferic. The L-value 

and values of n or N can be read directly from the position of the 
eq T 1 

whistler nose on the spectrogram. This method neglects corrections for 
dispersion in the conjugate ionospheres (Chapter 3) and for sub-ionospheric 
propagation t ime’ ( Chapter 4). This is not very serious, however, because 
the two corrections tend to cancel one another and because of the rela- 
tively poor precision inherent in this method. As mentioned in the pre- 
vious section, n^ depends sensitively on electron concentration models and 
woulc require a large number of overlays. It is therefore preferable to 

obtain n from n using Figure 12. 

1 eq 

An important advantage of this method is that it permits a quick 

identification of any structures in n or N profiles from a train of 

eq T 

nose whistlers. For example, if n is constant with L in some part of 

eq 

the magnetosphere penetrated by whistlers, the train of corresponding 

noses will follow an = constant curve in Figure 22. Electron concen- 

eq 

trat.on profiles in the magnetosphere frequently exhibits complex struc- 
ture!? such as those reported by Park and Carpenter [1970], These struc- 
ture!? can be easily identified on whistler spectrograms with the aid of 
the Dverlays. 

G . EMPIRICAL FORMULAS 

As shown in Figures 14 through 21, K, K , and K ^ vary only slowly 
with f T ^. These quasi-constants can be approximated by various curve 


51 


SEL-71-058 



fitting techniques to achieve desired accuracy. For slide rule calcula- 
tions, they can be approximated by constants. For the DE-1 model, the 
following approximate formulas can be used: 


Formulas 


f = 2.7 f f 
Heq n 


7 8.74x10 \ 
f Heq 


1/3 


f ’ t' 
n n 


n „ = 24 
eq \ T 5 


/f 't' 

N = 8 . 6x10 ( n — 
T \ L 


/f ' t’ 

”, ■ 


and for the CL model, 


Formulas 


f = 2.3 f T 
Heq n 


78.74x10 \ 
L = (— ) 


1/3 


Heq 


/f f t' 

, / n n 
n = 101 — 

eq V 5 


/f ’t 1 

» T - — w 


Maximum Error 
2.5 < L < 7 


1 . 1 % 


0.4 


7.6 


5.4 


25 


Maximum Error 
2.5 < L < 7 

7.6% 


2.5 


35.9 
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These- approximations by constant coefficients are reasonably good for 
diffusive equilibrium model, but they lead to somewhat larger errors in 
the case of the collisionless model. In particular, an attempt to use a 
constant coefficient for the collisionless model may lead to errors as 
large as a factor of 4, 

For machine processing, a parabolic approximation of the following 
form gives high accuracy with a reasonable amount of computation. 

K = a + a, F + a. F 2 (2.22) 

o 1 2 

where F = l°gio f n and ^ e atS are constants * The coefficients a 0 > a i’ a 2 
can oe obtained from a least squares fit to the values of the K’s in 
Tables 2 through 8. Tables 9 and 10 give the results of a least squares 
fit for the range 2.5 < L < 7 for the DE-1 model and R-4 model, respec- 
tively. The last column in the tables shows the maximum error due to this 
approximation in the range 2.5 < L < 7. 
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TABLE 9. COEFFICIENTS IN EQN. (2.22) FOR THE DE-1 MODEL 
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TABLE 10. COEFFICIENTS IN EQN . (2.22) FOR THE R-4 MODEL 
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III. CORRECTIONS FOR DISPERSION IN THE IONOSPHERE 


A. OUTLINE OF THE CHAPTER 

We use f and t to denote the whistler nose frequency and travel 
n n 

time at the nose observed after propagation through the ionosphere and the 
magnetosphere. Because the ionosphere has great variability, but rela- 
tively small effects on whistler propagation, we have chosen to account 

for the ionosphere by making ! corrections * to the observed f and t , 

n n 

obtaining quantities f f and t ' which would result from propagation through 

n n 

the magnetosphere only. Once f 1 and t T are obtained, the path latitude 

n n 

and the electron concentration along the path in the magnetosphere can be 

calculated by one of the methods described in the previous chapter. The 

purpose of this chapter is to examine in detail effects of the ionosphere 

on whistler propagation and to develop methods for obtaining f^ and t^ 

from f and t . 
n n 

B. WHISTLER DISPERSION IN THE IONOSPHERE 


The whistler travel time through the conjugate ionospheres, or the 
second term in Eq . (1.3) is 


t 


. (f) 

ci 


1 

2c 


/ 




ds 


where the subscript ci refers to the two conjugate ionospheres. The 
integration extends through both ionospheres below 1000 km altitude. Sev- 
eral simplifying assumptions and approximations will be used in this chap- 
ter in treating ionospheric effects on whistler propagation. Some of them 
are justified partly on grounds that the ionospheric contribution to 
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the total whistler travel time is small (typically less than 10%) . 
f 

Since *— « 1 in the ionosphere, we can write 
H 


D . 
ci 


t . f 
ci 


1/2 ~ L . f 

~ 2c J 

Cl 



(3.1) 


The quantity D is called dispersion, and Eq . (3.1) shows that the disper- 
sion in the ionosphere is independent of wave frequency. The dispersion 
is additive, so that 


D . = 
ci 


(D ) 

i northern 
hemisphere 


+ 


(°i > 


southern 

hemisphere 


In the remainder of this section, we will compute D from the ionospheric 
parameters . 

Assuming a horizontally stratified ionosphere and neglecting curva- 
tures of magnetic field lines, we write 


D . 

l 


1 

2c sin 6 




(3.2) 


where 6 is the magnetic dip angle, and h is the altitude. If we further 

neglect small changes in latitude as we move along inclined magnetic field 

lines, f can be expressed as 
H 


f 


H 


Ho 


+ 


h 

r 

o 


where f is the electron gy rof requency at h = 0 and at 
Ho 

interest. Substituting this into Eq . (3.2), we obtain 


the latitude of 
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1 


dh 


(3.3) 


D. 

1 


2c f 


T]2 

Ho 


sin § 



Equation(3 .3) has been evaluated numerically for several electron concen- 

0 

tration profiles at L = 4 where f = 1 . 57 v 10 Hz and sin6 = 0.957. The 

Ho 

ionosphere was assumed to be a simple cfChapman layer described by 

n = n max eXp | \ t 1 “ z " exp(-z)]j (3.4) 

h-h 

o 

z = — 

where n is the maximum electron concentration and h the height of the 
max o 

maximum concentration. Three different profiles corresponding to scale 

height H = 50, 75 and 100 km have been considered. Figure 24 illustrates 

the three profiles with h = 300 km and n adjusted so as to involve the 

o max 

same height-integrated content. The results are summarized in Table 11. 

Table 11 shows that the dispersion does not depend sensitively on the 

shape of the electron concentration profile. A comparison between the 

constant n case (A, B and C) and the constant integrated content case 
max 

(D, E and F) shows that the dispersion varies more with the content than 

it does with n . In estimating D it is therefore preferable to use 
max i 

integrated columnar content information obtained from Doppler shift mea- 
surements or Faraday rotation measurements of satellite beacon signals. 

Since about two thirds of the integrated content is above the F layer 
peak, topside sounder results are preferred over bottomside sounder results. 
The desired information on the ionosphere near a whistler path of interest 
is usually not available, particularly in the case of satellite data, and 
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1 
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ELECTRON CONCENTRATION (cm' 


TABLE 

11. IONOSPHERIC 

DISPERSION 


Scale 

n 

Columnar 

D. 

max 

Content 

l 

height 

(km) 

(el cm ) 

— 2 

(el cm ) 

1/2 
(sec ) 


^ -iO 

A 50 10 2.06 x 10 4.43 

B 75 10 6 3.08 x 10 13 6.31 

/? 1 q 

C 100 10 4.01 v 10 7.72 

D "50 10 6 2.06 x 10 13 4.43 


5 13 

E 75 6.68 x 10 2.06 x 10 5.17 

F 100 5.13 x 10 5 2.06 x 10 13 5.53 
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an ’average 1 ionosphere must often be used. 

The following formulas are offered for estimating from columnar 
content N or the F2 layer critical frequency f F2 . 


D. =1.15 N 1 / 2 
1 


(3.5) 


D. = 0.7 f F2 (3.6) 

l o 

12 2 

In the above equations, N is in units of 10 el/cm and f F2 in MHz. A 

7 o 

typical value of D ci (for two conjugate ionospheres) may be 8 during the 

day and 4 during the night near solar-cycle minimum. The corresponding 

values for solar cycle maximum may be 16 and 8. 

The values of f and sin6 in Eq. (3.3) vary with magnetic latitude. 

Ho 

However, in the principal range of whistler propagation 2 < L < 7 they can 
be regarded as constants, because their variations are not important com- 
pared to the uncertainties in electron concentration profile usualLy en- 
countered. The height of the F layer peak makes little difference in the 

calculated value of D. . A change in h by 200 km results in a change of 

i o 

only 5% in . 

C. CALCULATIONS AND RESULTS 

The whistler time delay is written as 

A 

t(f) = en 1 ^ 2 / 3 ( 0 , f ) At + D . f 1/12 (3.7) 

eq J 0 r ci 

The first term on the right hand side is the magnetospheric part from 
Eq . (2.19), and the second term the propagation delay in two conjugate 
ionospheres. Equation (3.7) is evaluated numerically by the same procedures 
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used in Section 2C in search of the minimum time delay. In the present 

case, however, f depends on n and D . as well as on the whistle? path 
n eq ci 

latitude. The dependence of f on n for constant L and D is illus- 

n eq ci 

trated schematically in Figure 25. The two dashed curves A and B repre- 
sent the first term in Eq . (3.7), or the travel time through the magneto- 
sphere. The two whistlers propagate on the same magnetic shell, and hence 

exhibit the same f 1 . The travel time for curve B is assumed to be 9 times 
n 

as large. The curve marked C represents the travel time through conjugate 

ionospheres, or the second term in Eq . (3.7). The solid curves are the 

result of adding "C to A and B. Evidently the effect of the ionosphere on 

f is less for larger magnetospheric electron concentrations. It is also 

apparent that the larger the value of D^, the larger its effect will be 

on f . 
n 

The calculations of f are repeated for variations in either D . , t \ 

n ci n 

or L while the other two are held constant. The values of D . used are 

ci 

D . = 2,4,8,16 and 32. The results of such calculations for L = 2.5, 4 

and 6 are shown in Figures 26 through 28 where the quantity (f -f^)/f^ is 

plotted against D . for several values of t T . The quantity (t -t f ) is also 
ci n n n 

calculated from Eq . (3.7) and plotted in Figures 29 and 30. 

D. EMPIRICAL FORMULAS 

Although f’ and t 1 can be obtained from f , t and D . with the aid 
n n n n ci 

of graphs illustrated in Figures 26 through 30, such a process is obviously 

cumbersome and time-consuming. In this section, we construct a few simple 

empirical formulas to approximate the results of the previous section. 

We expect from the earlier discussion that the difference between f^ 

and f f would increase with D . and decrease with t , so we try an equation 
n ci n 

of the form 
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with increasing electron concentration in the 



IODEL 

ODEL 



Figure 26. A plot showing the effects of ionospheric dispersion on the nos 
frequency of a whistler propagating along an L = 2.5 field line 
for different values of t' and for two different field-line dis 
tribution models. 
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Figure 27. A plot similar to Figure 26 but for an L = 4 field line 




0.6r~ L=6.0 
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Figure 28. A plot similar to Figure 26 but f 
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igure 29 but for 




f ’ 
n 


f - 


- y' 


D . 
ci 

t 


n 

It is found that y f is 

decreases with increasi 

tional to f , we write 
n 


n 

nearly constant for a given 
-1/3 

ng L. Since L or f , and 
Heq 


L-value, but that it 

f T is roughly propor- 
Heq 


f - f f 
n n 


f 1 


n 


Y 


D . 

Cl 



(3.8) 


Equation (3.8) with a constant y is a good approximation to the results 

of the previous section for a wide range of f , t and D . Empirically, 

n n ci 

the best value of y is found to be 0.17 for DE models and 0.15 for the CL 
or R-4 model. Thus, 


f 1 = 


n 1 + 0.17 (D It f Y 3 ) 
ci' n n 


(3.9) 


for DE models, and 


1 + 0.15 (D 


,/t f 'Z 3 ) 

iv n n 


(3.10) 


for the collisionless model. The error in due to the above approxima- 
tion is plotted in Figures 31 and 32 against D for a range of values of 

t r a nd L. The fact that the per cent error is generally larger for larger 
n 

D and for smaller t 1 means that the approximation by Eqs . (3.9) and (3.10) 
ci n 

is better when the ionospheric part of the whistler dispersion is relatively 


small.. For all realistic values of D . and t\ the error in f f is less 

ci n n 

than ^ 3%. Several curves that surpass the 3% error level in Figures 31 
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DE-1 MODEL 



70 


Figure 31. A plot of the percentage error in 
for a range of values of L and t T 




71 


of values of 


and 32 correspond to extremely low electron concentrations in the magneto- 


sphere. For example, the solid curves marked t' - 0.2 sec are approaching 

n 

the theoretical low limit of 0.18 sec which is the propagation time along 
the length of an L - 4 field line in free space. When electron concentra- 
tions are that low, errors due to other sources become more important than 
the error due to the empirical formulas given by Eqs. (3.9) and (3.10) (see 
Chapter 6) . 

We now consider expressions for t f . From Eq . (3.7), 

n 


t t 1 + D . f 
n n ci n 


'V 2 ... + D , f 1 2 

n ci n 


Figure 33 is similar to Figure 25 except that the ionospheric dispersion 


is exaggerated in order to illustrate details near the nose frequency. It 

, - 1 ' 2 

is evident in the figure that t - t f is slightly less than D . f* 

n n ci n 


but slightly greater than f 


- 1/2 


A compromise is thus made: 


, /f + f’ \ -1/2 

t 1 = t - D. / n n] 

n n ci ^ 5 / 


Equation (3.11) approximates the results of the previous section with less 

than 1 msec error for wide ranges of f , t and D . . 

n n ci 
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Figure 33. A sketch showing the effect of ionospheric dispe 
the minimum time delay of a nose whistler. 









IV. SUBIONOS PHERIC PROPAGATION 


The whistler travel time is usually measured with respect to the 
causative sferic on a whistler spectrogram, and it is necessary to account 
for the travel time of the signals in the earth-ionosphere waveguide. 

Figure 34 illustrates the geometry and the appearance of the sferic on a 
frequency-time spectrogram. At the top of the figure, R and T represent 
the locations of the whistler receiver and the sferic, respectively. A 
sferic often shows a small dispersion and a sharp low frequency cutoff as 
illustrated at the bottom of the figure. The cutoff occurs at the earth- 
ionosphere waveguide cutoff frequency which is 1.67 kHz for a lossless 
waveguide with a height of 90 km. The reader is referred to Helliwell 
[1965] for more detailed discussion of the dispersive effects of the 
earth-ionosphere waveguide. In this chapter, we estimate the subi onospheri c 
propagation time of the lightning signals that reach the whistler receiver. 
We will use t to denote the whistler travel time through the ionosphere 
and the magnetosphere and t to denote the travel time of the whistler mea- 
sured on a spectrogram with respect to the sferic. 

We consider a simple case in which the lightning, whistler duct and 
the receiver are all in a magnetic meridional plane. The distance traveled 
by the direct wave can be written 

rtr 

l = — 2- (0 + 0 ) (4.1) 

D 180 'T - R 

where 6 and 6 are the magnetic latitude in degrees of the sferic and the 
k T R 

receiver, respectively. The sign of the latitude is positive in both \he 
northern and the southern hemispheres. The subi onospheric distance tra- 
veled by the whistler wave is 
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(4.2) 


irr 

= 180 ( ~ ^dI + I^R “ ^ 

w here j) is the latitude of the end points of the whistler duct. The 
correction necessary for the subionospheric propagation is then 


^W 


T - 


where c is the speed of light in a vacuum, 
the following equations are obtained: 


(4.3) 

Using Eqs . (4.1) and (4.2), 


t 

- j — 6 .65 

X 

10 4 

K 

sec 


for 

0 

> 

K 

a nd 0 > 

R — 

0 

(4.4) 



-4 






r D 




and 0 < 0 


t 

- T = 6.6 5 

X 

10 

K 

sec 


for 

0 T 

> 

K 


(4.5) 


- T = 6.6 5 


-4 











t 

X 

10 

K 

sec 


f or 

0 


K 

and 0„ > C 

S 

(4.6) 







r T 



D 

- 

T= 6.65 X 10" 

X 

+ 1 

v- 

0 ) 

for 

0 

X \ 

0 

and 0 R < ? 

i 

D 

(4.7) 


For whistler receivers at high latitudes, 0^ is usually larger than 
0 , and therefore, Eq . (4.4) or Eq . (4.6) would apply. For 0^ = 60° 

(L = 4) , t - t may vary from 0(0 = 0) to 0.04 sec (0 = 60°) depending on 
0^ Since 0^ is usually not known, a constant value of t - 0.03 sec 

may be used for all medium to high latitude whistlers as a first order 
correction . 

If we remove the assumption that the sferic and the duct lie in the 
same meridian plane as the receiver, the distances i and l should be 
great circle distances. However, in view of the fact that no means of 
identifying the location of sferics and the longitude of whistler ducts 
are available at present, we will not consider such refinements here. 
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Figure 34. The sketch at the top shows a whistler mode and a wave- 
guide mode propagation paths of lightning signal from 
its source T to the receiver R in the conjugate hemisphere. The bot- 
tom sketch shows how the received signal appears on a f requency- t ime 
spectrogram. The wave-guide mode signal shows the effect of the wave 
guide cutoff at low frequencies . 



V. SUMMARY 


In this chapter, we briefly summarize the procedures for determining 

the path latitude and electron concentration along the path from the 

wh:.s:ler parameters f and t . 

n n 

1. Scale f and j from whistler spectrogram. 

n n 

2. Make corrections for subi onospheri c propagation as discussed in 

Chapter 4 to obtain t . Typically, t = t + 0.03 sec. 

n n n 

3. Decide on a magnetospheric electron concentration model to be 
used . 

4. Estimate dispersion in the two conjugate ionospheres D . as dis- 

ci 

cussed in Chapter 3. 

5. Make corrections for the ionospheric dispersion to obtain f^ and 

t f as follows: 
n 

(a) for diffusive equilibrium models (DE-1 through DE-4) 


f 1 
n 


1 + 


0.17 D A 

Cll 


t f 
n n 


173 


(5.1) 


f + f ' \ ^ 

n n 


t T = t - D . 
n n ci 


(b) for collisionless models (CL and R-4 models) 

f 


n 1+1 


0.15 D .n 
ci 


t f 
n n 


173 


t T = t - D . 
n n ci 


-l/2 



+ f 1 
n 

2 


(5.2) 


(5 3) 
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6. L, n , N and n can be obtained from f' and t’ bv one of the 
eq T 1 n n 

following methods depending on the desired precision and availa- 
bility of a computer or a calculator. 

(a) L, n and are read off the graphs in Figures 7, 8 and 9 

for DE models and Figures 7, 10 and 11 for the CL model or R-4 

model. n is obtained from n using Figure 12. 

1 eq 

(b) Calculate L, n , and us: *- n £ the following formulas 


f = K f f (5.4) 

Heq n 



(5.5) 


(5.6) 


(5.7) 


(5.8) 


using the values of K, K , and obtained from Figures 14 through 21. 

(c) Use Eqs . (5.4) through (5.8) above, but obtain the values of 

2 x 

K 1 s from empirical formulas of the form K - a^ + a^F + a^l whert' 

F = log f’. The coefficient a's are given in Tables 9 ami 10 for 
10 n 

the DE-1 model and R-4 model . 
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VI. UNCERTAINTIES IN THE WHISTLER METHOD 


A. INTRODUCTION 

In this chapter, we discuss uncertainties in the whistler method des 
cr.ibed in this report. A few general remarks will be followed by more de 
tailed examination of various sources of error. 

As mentioned earlier in Chapter 2, whistlers are most sensitive to 
the conditions at the top of the propagation path. To illustrate this, 
we return to Eq. (1.2). At the nose frequency, 


i r p 

r 175 J ./2„ V 3 ' 2 

" J V (1 - - > 


ds 


From the results of Chapter 2, f ^ 0.377 f for diffusive equilibrium 

n Heq 

models, so t is approximately proportional to 
n 


/ 


~~ 7 f 2 vT 3/2 

f 1/2 (1 - 0.377 

H I 

H 


ds 


Figure 35 is a plot of the value of the integrand, 

f 


I 


, 1 f \ 3 / 2 

r H 1/2 1 - 0.377 -fSS) 


V 

normalized to the equatorial value, as a function of distance along lines 
of lorce for several L- values. The DE-1 Model was used for the electron 
distribution along field lines. The terminal point of each curve corre- 
sponds to 1000 km altitude, and the circle on each curve to 30° dipole 
latitude. The whistler time delay is proportional to the area under the 
curve, and it can be seen that nearly 80% of the time delay occurs within 
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— 30° of the equator. The sharp f hooks f in the curves are due to rapid 

increases in electron concentration near the base of the magnetosphere. 

If a diffusive equilibrium model with different parameters (temperature 

and ionic composition) is used, there are changes in the shape of the 

’hocks’ , but little effect at distances greater than 1000 km from the base 

of the magnetosphere. For collisionless models, f a? 0.41 f , and the 

n Heq 

has been calculated at several 

L--values and is plotted in Figure 36. In this case, — 60% of whistler 
propagation delay occurs within 30° of the equator. Because of these 
facts, electron concentration near the equator as measured by the whistler 
technique is remarkably independent of assumed electron distribution along 
field lines. On the other hand, electron concentration near the base of 
the magnetosphere depends sensitively on assumed models. 

Information on electron concentration in the magnetosphere usually 
involves extrapolation along magnetic field lines by assumed models. For 
example, electron concentration measured near the base of the magnetosphere 
by incoherent backscatter radars, satellite-borne sounders or probes can 
be extrapolated to the equatorial plane. When computing magnetospheric 
tub? content, however, the whistler technique has the advantage of measuring 
electron concentration accurately near the equatorial plane where most of 
the tube volume lies. 

We now examine more quantitatively various sources of uncertainties 
in the whistler method. 

B. UNCERTAINTIES DUE TO THE ASSUMED ELECTRON CONCENTRATION MODEL 

For given f f and t ’ , Eqs. (5,4) through (5.8) can be written in dif- 
n n 


f 
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N 
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Figure 36. A plot similar to Figure 35 but for the CL model of field-li 
distribution of plasma. 



ferential form as follows: 


m 


^f 

Heq 



The 1 s in the above equations represent differences between electron 

concentration models f and their values can be obtained from Figure 1 A 

through 21. For a whistler with observed f , f’ depends somewhat on the 

n n 

electron concentration model used (see Eqs . 5.1 and 5.3), but this is a 

sma'.l effect and will be ignored. An example given below illustrates the 

effects of changing the model from DE-2 (3200 °K, 90% 0 + , 8% H + , 2% He + ) 

to DE -3 (1600 °K, 50% 0 + , 40% H , 10% He + ) for a whistler with f 1 - 5 kHz. 

n 

From Figures 14 through 17, we obtain 


AK 

K 


0 . 012 , 


AK 


eq 


eq 


0.071, 


AK 


A K, 


0.0 27 and 


0 17 


83 


SEL-71-058 



Substitution of these values into the above equations yields 


AL 

L 


0 .004, 


An 

— — = - 0.091, 
n 

eq 


AN 

T 


0.031 



0.19 


It is seen that N and n are relatively insensitive to changes in model, 
l eq 

whereas n^ is more model-dependent. 

C. MEASUREMENT ERROR 

We assume that measurement errors in f and t result in the same 

n n 

percentage errors in f 1 and t f . Since the K’s in Eqs . (5.4) through (5.8) 

n n 

are only very slow functions of f^, we further assume that the K f s are 

constant within the range of measurement errors in f . Equations (5.4) 

n 

through (5.8) can then be written in differential form as 


Af 


Heq 

Heq 


* f ; 


f 1 
n 


AL 

L 


1 

3 



( 6 . 1 ) 


A n 


eq 


n 

eq 





( 6 . 2 ) 



(6.3) 


(6.4) 
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In routine scaling, measurement errors can be easily maintained below 


^ 3% for f T and ~ 1% for t 1 . Equations (6.1) through (6.4) show that a 
n n 

1 % error in t* causes a 2% error in computed n , N and n . An error of 
n eq T 1 

3% in f f causes a 1% error in the computed L-value, a 4% error in N and 
n T 

8% errors in n and n . With high quality whistler traces and sufficient 
eq I 

care in scaling, measurement accuracies can be improved over the figures 
quoted above, 

D. UNCERTAINTIES DUE TO THE IONOSPHERE 

Uncertainties in ionospheric dispersion appear as uncertainties in 

f' and t'. From Eqs . (5.1) and (5.3), we can write 
n n 


A* f 

n 


-0.17 


t f 
n n 


173 


AD, 


Cl 


for DE models and 


f\f 1 
n 


f f 


n 


-0.15 


t f 
n 


V 3 


AD . 


Cl 


for collisionless models. 


From Eq . (5.2) 


At’ 
^ n 


t 1 
n 


1 


t f 
n 


T72 


AD 


Cl 


for both DE models and 
collisionless models. 


Equations (6.1) through (6.4) 

in L, n , N m and n . 
eq T 1 

As an example, consider 

(DE model) with f = 6 kHz (L 
n 

u -1/2 

D by 1 sec > results in * 
ci 

these values into Eqs. (6.1) 


can then be used to estimate uncertainties 


a whistler propagating in 

~ 4) and t = 1 sec. An 
Af’ n At; 

— - = -0.009 and = - 

f t 

n n 

through (6.4), we obtain 


the plasmasphere 
overestimate of 
0.002. Subs ti tut ing 


85 


SEL-71-058 



AL 

L 


0 .003, 





In addition to uncertainties in D ci , the empirical approximation of 

Eqs . (5.1) and (5.3) introduces small errors in f f . The magnitude of this 

n 

error depends on magnetospheric parameters and D . (see Section 3D) and can 

ci 

be estimated with the aid of Figures 31 and 32. 

E. UNCERTAINTIES DUE TO SUBIONOSPHERIC PROPAGATION 

Whistler propagation in the earth-ionosphere wave guide has negligible 

effects on f^ unless f^ approaches the wave guide cutoff frequency of 

^ 1.6 kHz. Corrections for subionospheric propagation, however, must be 

made to time delays measured on a whistler spectrogram. As discussed in 

Chapter 4, the amount of correction necessary is typically ~ 0.03 sec with 

uncertainties of ~ 0.015 sec due to uncertainties in the location of light- 
en 

ning. For t =1.5 sec, this corresponds to — r ~ 0.01, and from Eqs. 
n t 

An AN An. n 

(6.2) through (6.4) , = — - = -0.02. 

n N n„ 

eq T 1 

F. SUGGESTIONS FOR FURTHER WORK 

A number of simplifying assumptions were made in this study in order 
to facilitate calculations and arrive at simple recipes for routine process- 
ing of whistler data. It was also necessary to use several models for 
electron distribution along geomagnetic lines of force. When a need arises 
for improved accuracy in whistler methods, these models and assumptions 
must be re-examined critically. 

As mentioned previously, there is empirical evidence in support of a 
diffusive equilibrium model inside the plasmapause under normal conditions. 
Outside the plasmapause, however, the idealized collisionless model used 
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here is probably extreme, and the actual situation may be better repre- 
sented by a model falling between the collisionless model and the diffu- 
sive equilibrium model (see Section 2B) . Banks, et al [1971] recently 
pointed out a theoretical possibility that a low plasma pressure outside 
the plasmapause following a magnetic storm may cause the ionospheric plasma 
to flow upward at supersonic speeds. The field line distribution of plasma 
may be very complicated under such circumstances. More detailed theoreti- 
cal and experimental studies are needed to improve the electron concen- 
trat ton models used outside the plasmapause. 

The expression for refractive index used in this study is for a cold 
plasma., and it is further simplified by neglecting the additive 1 i n Eq . 
(1.1) on the assumption that plasma frequency is much higher than electron 
gy rof requency . These assumptions need further critical examination outside 
the plasmapause, where electron concentrations are low and the temperatures 
are high . 

Another important simplification made in this study is the use of a 
dipcle magnetic field. At large L values and during magnetic disturbances, 
the earth’s magnetic field is distorted significantly by the solar wind, 
and the ring current. A detailed investigation should be made of whistler 
propagation in distorted magnetic fields and the possibility of detecting 
such distortions by the whistler technique. 

An improvement is desired in estimating subionospheri c propagation 
time and correcting for the effect of wave guide dispersion on nose fre- 
quency when the nose frequency is close to the wave guide cutoff frequency. 


87 


SEL-71 -058 



APPENDIX - DIPOLE GEOMETRY 


A. MAGNETIC FIELD STRENGTH 

In a centered dipole approximation, the earth's magnetic field 
strength is given by 

B = 0.312 (r^/r) 3 (l + 3 sin 2 ^) 1 ^ 2 gauss (A. I) 

where 6 is the latitude, r the geocentric distance and r the mean radius 
r o 

of the earth, 6370 km. Figure 3 illustrates the dipole geometry and the 
symbols used. The electron gyrof requency is given by 


f 


H 


1_ Be 

2tt m 


5 3 9 1 / 2 

8.736 X 10 (r /r) (1+3 sin 0) Hertz 

o 


(A. 2) 


Figure 37 shows the equatorial field strength B (and f ) 

eq Heq 

r , while Figure 38 shows the latitudinal variation of the 
eq 

at the earth's surface B . 

o 


plotted against 
field strength 


B. L- VALUE 

A dipole field line is described by 


2 , 
cos 0 


2 , 

cos 0 


(A. 3) 


where the subscript o refers to the earth's surface. 
Mcllwain’s L parameter is 


L = 



o 


2 , 

cos 0 


o 


In a dipole, 


(A. 4) 
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Figure 37. A plot of the earth’s magnetic field strength and electron 
gy rof requency at the equator as a function of L according 
to a centered dipole model . 
























If the reference level is changed to 1000 km altitude, 0^ can be substituted 
for 0 o by using Eq . (A. 3) and we obtain 

r 

1 1 

L = — (A * 5) 

o cos 0 

Figure 39 shows plots of 0 and 0 versus L-value. 

o 1 

C. DIP ANGLE 

The radial and tangential components of B are given by 


B 

r 


0.624 ( r /r) cos( 
o' 


B . = 0.312 

e 



sing 


These vector components are related to dip angle by 


(A. 6) 


B 

tang = — = 2 cot9 = 2tan0 (A. 7) 

B e 

Figure 40 is a plot of 6 as a function of 0. The horizontal scale at top 
is the L-value at the earth’s surface. 


D. LENGTH OF FIELD LINES 

The arc length along a dipole field line between the equator and 6 
can be written as (see Helliwell [1965]) 



r 

o 


cos 



(x + sinhx coshx) 


(A. 8) 


where sinhx = /3~ sin0. By substituting sinhx 


(e X -e X ) / 2 and 


coshx = 
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Figure 40. A plot of the dip angle as a function of latitude for a dipole field 
horizontal scale at the top shows dipole L parameter. 



X ~ - X 2 9 

(e + e ) / 2 , and using the identity cosh x - sinh~x = 1. Eq . (A. 8) can 
be rewritten as 


1 

2j2 


2 / 
cos 0 


’n( /T 


sin0 + / 1 + 3 sin 




4/3 sin0 '1 


4 


2 

s i n“0 


(A. 9) 

Figure 41 shows a plot of the arc length measured from the dipole equator 
as a function of latitude for several L-values. Figure 42 shows the arc 
length from the equator to the 1000 km level plotted against L. 

E . TUBE VOLUME 

2 

We define V as the volume of a tube of force with 1 cm cross- 
sectional area at 1000 km and extending to the equator. We write 

equator 

V T= J |i,ls 

1000 km 

and use Eqs. (A.l) and (2.15) to obtain 



o 


The integration yields 
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Figure 41. A plot of the arc length of the earth's dipole field lines as a 
function of latitude for several L values. 
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V = 


4 4 

r L 
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(1+3 sin 2 <f)) 1 ' 2 


3 5 7 7 

64 Sin ^l + 64 Sin3 ^l + 320 
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448 


sin70 


Figure 43 shows a plot of 


versus L. 
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